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Binary black hole coalescence has its peak of gravitational wave generation during the “plunge,” 
the transition from quasicircular early motion to late quasinormal ringing. Although advances in 
numerical relativity have provided plunge waveforms, there is still no intuitive or phenomenological 
understanding of plunge comparable to that of the early and late stages. Here we make progress 
in developing such understanding by relying on insights of the linear mathematics of the particle 
perturbation model for the extreme mass limit. Our analysis, based on the Fourier domain Green 
function, and a simple initial model, point to the crucial role played by the kinematics near the 
“light ring” (the circular photon orbit) in determining the plunge radiation and the excitation of 
QNR. That insight is then shown to successfully explain results obtained for particle motion in a 
Schwarzschild background. 


I. OVERVIEW 
A. Background 

Interest in the coalescence of binary compact objects was motivated by the development of gravitational wave (GW) 
detectors, and the need for a better understanding of the endpoint of binary merger [1]. Especially in the case that 
such mergers end in black hole formation these might be the most powerful, and certainly the most interesting sources 
of detectible GWs. Beyond that main motivation, there are several other reasons such mergers are of interest to 
the physics research community; see Ref. [2] for a detailed review. From the outset it was understood[5] that binary 
coalescence, driven by GW radiation, could be viewed as having three stages. The first is the gradual inspiral of the 
objects in orbits, most probably nearly circular orbits, gradually decreasing in orbital energy, and hence in radius, 
driven by the radiation reaction to the GWs being emitted. This first stage can be understood as a quasi-Newtonian 
process with GW radiation reaction treated as a time averaged loss of orbital energy. Considerable accuracy could be 
achieved by using post-Newtonian (slow motion, weak field) approximations [4]. 

The last stage of binary coalescence is the oscillation of the black hole formed by the coalescence. These oscillations, 
after initial transients, are the “quasinormal ringing” (QNR) of the final hole. It has been known since the early 
1970s[S] that these QNRs constitute a spectrum of complex (i.e., damped) frequencies characteristic of the mass and 
spin of the final black hole. 

The second stage, the middle stage in the binary coalescence, may usefully be called the “plunge.” It is the end of 
the early gradual inspiral, and determines the magnitude and phase of the QNR. Most important, this relatively short 
transition is the stage during which the major part of the GW energy is radiatedjS]. Ironically, this most-important 
stage is the part of the coalescence that is least understood. 

For a wide variety of merger scenarios, the excitation of QNR is known. The remarkable progress in numerical 
relativity in the past decade[7] has given detailed results for the GW signals produced. These results, however, do 
not give a qualitative, phenomenological understanding of the plunge process, an understanding at the same level as 
that we have for the early inspiral and late QNR. In addition to the satisfaction that such an understanding could 
provide, there are possible practical advantages. With an intuition about the plunge, computational searches could 
be guided to particularly interesting scenarios. In this paper we attempt to give a first step, an intersting first step, 
towards the understanding of the plunge. 


B. Introduction and Summary 

We attempt here to give a simple understanding, more precisely a simple view of the generation of plunge radiation 
and the excitation of QNR. The ultimate goal is a set of guidelines that could allow a complete quasi-analytic 
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description of merger. It should be noted in this connection that the difficulty posed by the plunge, the end of inspiral 
and beginning of QNR, has been recognized as a challenge in the effective one body formalismjS]. 

The result we present is a very simple general picture: the plunge is determined by the conditions of the binary 
associated with the “light ring,” (LR) the location of photon orbits. This picture is supported both by a mathematical 
reason that the picture should apply, and by a variety of examples. The simple general picture, however, must be 
viewed as somewhat speculative since it depends on several approximations. By far, the most important of these 
is the particle perturbation approximation. We have also limited our analysis to nonrotating (Schwarzschild) black 
holes and, in some aspects of the mathematical analysis, to simplified models. Despite these approximations, we feel 
that the picture that emerges is very compelling as a first step that can motivate and guide further work. 

There has been considerable past work[9UI^ reported on the connection between the LR and QN frequencies. Of 
particular relevance to our current work is the analysis of the multipole content of merger ringdown radiation and 
the relationship of the start of ringdown to the passage of the particle through the light ring[TS]. Our interest here 
is somewhat different, not only to observe the connection in the radiative signal but to understand its origin, and to 
understand the transition to the ringdown and the peak of the radiation. 

The approach presented in the sections below consists of the following steps. First, in Sec. |Tl]we introduce the 
particle perturbation approximation and the idea of using the Fourier domain Green function (FDGF) to understand, 
with simple mathematics, the analytic connection between the motion of our point particle and the excitation of 
plunge radiation and QNR. With this technique, the radiation is reduced to a single, relatively simple, integral over 
the path of the particle. The use of the FDGF, in fact, is one of the elements of our hope that the plunge can be 
understood. The basis of that hope is that the QN frequency is a pole in the FDGF. This is true in very broad 
generality, so the way in which a particle source “excites” a pole should have significant universality. Excitation in a 
simple model may therefore tell us a great deal about a much more complicated realistic model. With this motivation 
we introduce, in this section, a model with a very simple FDGF. Even this simple model, however, is complicated 
enough in its details to be distracting, so most of these details are deferred to an appendix. 

Next, in Sec. Ill we define a class of trajectories that allow for simple analysis combined with sufficient flexibility 
to test models. In this section, also, we show that the perturbations computed from the FDGF, for the simple model, 
agree to high accuracy with numerical evolution of initial data. We then, in Sec. |IVj use the FDGF to investigate 
what it is about the trajectory that determines the excitation of QNR. Our investigations in this section are limited to 
radial motions, since our model is based on a scalar charged particle source of a scalar field. For the scalar case, orbital 
motion modifies the source with a multiplicative, time dependent angular factor. By contrast, for the gravitational 
and electromagnetic case, the source term has extra additive terms for angular motion. Our scalar model, therefore, 
will not give results that are representative of gravitational wave generation, it does however give interesting insights 
into the excitation of QNR and to the plunge radiation. 

Our conclusions from Sec. [IV] are extended and tested in Sec. |Vj where we show them to be valid for trajectories 
in the Schwarzschild background. In this section we include investigations both of radial motion and motion with 
angular orbital motion and we show in this section that the plunge radiation and QNR are determined by the angular 
as well as the radial motion of infall. 

In Sec. lYll we give a viewpoint on why the LR plays the important role in the generation of QNR. We conclude in 
Sec. VII and relate the work presented to directions for further investigations. Some details of the complete FDGF 
are given in Appendix [A| 

We use units in which c = G = 1 and other conventions of the text by Misner, Thorne and Wheeler [TB]. The 
following acronyms will be used throughout the paper: QN=quasinormal; QNR=quasinormal ringing; LR=light ring; 
FDGF=Fourier domain Green function. We will typically use the capital T to indicate particle time, and the lower 
case t for other uses, in particular as part of the retarded time of received radiation. 


II. PARTICLE PERTURBATION APPROXIMATION AND FOURIER DOMAIN GREEN FUNCTION 

A. Particle perturbation model 

We start by considering a model of a perturbation field $ generated by particle motion in a Schwarzschild back¬ 
ground. In effect, this particle perturbation model is justified by the fact that it applies in the limit of an extreme 
ratio of the masses of the binary compact objects [T7ll20j . But an important question is whether phenomena in the 
extreme mass limit are qualitatively similar to the phenomena for comparable mass binaries. Many examples show 
that this does seem to be the case. One such example is the analysis of the antikick phenomenon found with numerical 
relativity computations for comparable mass |21lEB] . With particle perturbation analysis it has been shown 

that this paradoxical phenomenon has a simple underlying explanation. The current paper uses the history of such 
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successes as the reason to look for particle perturbation modelling as an appropriate step toward understanding the 
plunge and the excitation of QNR. 

The great advantage of particle perturbation theory is, of course, that the equations describing the fields are linear, 
and allow the insights of a Fourier transform and, for a spherically symmetric background, of multipole decomposition. 

For greatest simplicity we take our model to be based on a scalar perturbation field, <&, coupled to a scalar-charged 
particle. We start the analysis by decomposing d) into spherical harmonics r~^ ^tmYim and by writing the 
equation for the tm multipoles as 




dt"^ 




Stra{r*,t) . 


( 1 ) 


Here r* is the Regge-Wheeler [23 “tortoise coordinate” defined in terms of the standard Schwarzschild areal coordinate 
r, by 


dr/dr* = 1 - 2M/r . (2) 

The crucial feature of the r* coordinate is that it remaps the semi-infinite scope (2M to oo) of the Schwarzschild areal 
coordinate to — oo to oo. 

The “curvature potential” [23 , V^{r*) depends on the type of field (scalar, electromagnetic, gravitational 
perturbation) [27] . What is common to all these potentials is that for r* ^ M the potential has the form of a 
centrifugal potential. 


for r* > M , (3) 

and that as r* —> —oo the potentials fall off exponentially in r*, that is, as (1 — 2M/r). The transition between the 
flat spacetime centrifugal potential and the dramatically decreasing exponential form occurs around the “light ring,” 
r = 3M, the location of circular photon orbits. 

For the problem in which we are interested, the source term Sirn{f*T't) is a point particle moving in the equatorial 
{9 = 7r/2) plane with radial position r* = F{t). The nature of the source term, like the specifics of the potential, 
depends on details of the model. The source term may contain a delta function, or derivatives of delta functions, or 
sums of such terms. Our paradigmatic simplest case will be the scalar problem, with a source term that contains only 
a delta function, not its derivatives. In this case the source will have the form 


Stm{r*,t) = f{t)5{r* -Fit)) . 


( 4 ) 


Here /(t), which in general may depend on ^,m, can be used to represent, e.g., time dependence of scalar charge or, 
as shall be discussed later, orbital motion. This factor will almost always be taken as / = 1, but will be useful in 
Sec. HI for comparing numerical results. 

We now change notation in two ways. First, we drop the £,m indices which will be clear from the context of the 
models to follow. Second, we will want to generalize the background rather than be limited to that of the Schwarzschild 
geometry. Our generalization, in fact, focuses on Eq. Q as a wave equation in one spatial dimension. To emphasize 
this, we replace r* by the coordinate x, with the understanding that it ranges from —oo to -l-oo. Equation Q, then 
is replaced by 


^2^ 

9x2 


^2^ 


V{x)'i> = f{t)6 (x - F{t)) . 


( 5 ) 


B. The Fourier Domain Green Function 

We define the time domain Green function G(x, a]t — T) by the differential equation 

a2(-i 

^ ^ = Six - a)Sit - T) , (6) 

and the same boundary conditions that apply to Eq. Q: outgoing radiation (i.e., propagating to larger x) at x —> oo, 
and ingoing radiation (i.e., propagating to lower x) at x —> —oo. Since 

/ OO pOO 

/ 6ia — F[T])Six — a)6it — T) dadT, 

-OO j —OO 


( 7 ) 
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it follows that the solution to Eq. ([^ is 

/ oo noo 

/ G{x,a-,t-T)f{T)S{a-F[T])dadT. 

-OO ^ —OO 

Next we introduce the Fourier domain Green function (EDGE) Q by 

G{x,a;t — T) = — / e~''‘^^*~'^^Q{x,a;uj)duj, 

J — OO 

where (from Eq. Q), ^ satisfies 

^ + {uj'^ -V{x)) g = 5{x - a). 

We can put together Eqs. (|^ and ([^ to get 

-I pOO pOO pOO 

q,(t,x) = — / / e-^‘^^*-^^g{x,a;uj)f{T)S{a-F[T])dujdadT. 

J — OO —oo J — oo 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


C. The FDGF for a simple model 


The FDGF for the Schwarzschild curvature potential has a complicated analytic structure. The infinite set of 
quasinormal frequencies must appear, and - even worse - there is a branch point at w = 0 associated with the 
late-time power-law tails |28j. This complexity is a disadvantage for intuitive insights, and adds significant difficulty 
to evaluations of Green function solutions. For that reason we illustrate some of the features of our analysis with a 
model that has a very simple Green function. The truncated dipole potential (TDP)[5S] is defined by 


y — [ + 1)/^^ = 2lx^ for x > xq 

( 0 for X < Xq 


( 12 ) 


In this simple model, the sharp edge a.t x = Xq plays the role of the steep dropoff of the Schwarzschild potential for 
r less than around 3M. With the substitution of Eq. (12) for the Schwarzschild curvature potential, the solutions to 


Eq. (101 take an elementary form both for x < xq and x > xq- Matching conditions then lead to a simple polynomial 


equation to be solved for the QN frequencies. 

There are disadvantages to making the substitution of the TDP for the Schwarzschild potential, both obvious 
and not. The obvious disadvantage is that the substitution raises the question whether the insights provided will be 
applicable to the Schwarzschild background. A partial answer to this question is the observation that the phenomenon 
of complex frequencies occurs in the TDP (and in other truncated multipole potentials). It is reasonable to infer that 
the excitation of the TDP and the Schwarzschild potential will show qualitatively similar dependences on particle 
motions. The real answer to this question, however, is found in the models of Sec. |Vj the insights from the TDP 
model do seem to give correct insights into the nature of plunge radiation and QNR in Schwarzschild models. 

The not-so-obvious disadvantages of the TDP model are, however, worth a consideration. One of these disadvantages 
is the “quality factor” of the QN modes, the ratio of the real (oscillatory) part of a QN frequency to the imaginary 
(damping) part. For the least damped of the Schwarzschild quadrupole modes, the modes of prime importance to 
GW astrophysics, this ratio is 4.2, while for the single TDP pair of modes it is 1.0. Due to the low quality factor, 
QNR in TDP radiation does not have the same lightly damped sine wave appearance as for QNR in Schwarzschild 
models. Worse than this visual inconvenience is the fact that the damping part of the QNR for TDP models is smaller 
than that for the Schwarzschild models. (More specifically, the TDP damping frequency is l/(2a:o) or 1/6M if we 
associate xq with 3M. The Schwarzschild damping frequency 0.088965/M is about half of that.) We will argue that 
QNR is determined by the features of the orbit within a damping time of crossing the LR. The longer damping time 
for TDP than for Schwarzschild, means that this principle of LR importance is tested more stringently in the TDP 
models than in Schwarzschild models. 

It would be ideal, of course, to have a model for the potential term that has all the features we would like, but for 
which the FDGF has an elementary form. A sharp cutoff at some Xq, as in Eq. (12), guarantees that for x < Xq the 
Green function has a simple elementary form. The challenge is to find the right form of V(x) for x > xg. One such 
possibility is the Poschl-Teller potential [TUI, sech^x, but this potential does not approach the flat spacetime centrifugal 
potential at large x and hence is a questionable subsitution. An obvious question is why we do not use a higher order 
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truncated multipole potential, i.e., £(£ + l)/x^, with an integer i larger than unity. For all such potentials the Green 
function has only elementary functions. As £ increases, however, the Green function - though elementary- becomes 
more and more complicated. Furthermore, the quality factor of the least damped QN mode increases slowly with 
increasing 1. We would have to go to £ of order 10 to get a ratio of damping to oscillation that is qualitatively similar 
to QNR in the Schwarzchild geometry [ 30 ]. 

The complete FDGF for TDP, and its derivation, are given in Appendix [A} The explicit expressions for the FDGF 
depend on the relative location of a;, Xg, and particle position a, but all expressions are simple. As an example, for 
a < xq < X the expression is 


G 

where uji ^2 are the TDP QN frequencies 


2 (uj — uji){uj — uj2) 



2x0^ 2x0 ’ 


(13) 


(14) 


which show up, as expected, as poles in the FDGF. 

For a known trajectory x = F(T), we ca n find the solution to Eq. (§, for X > xo, at all times, by using G of Eq. ( [T^ 
for a < Xq, and the Green function (in Eq. (A21) for a > Xq. To do this we introduce notation for times of importance 
Fcross and Tq, important times along the trajectory. (The notation Tg arises from the step-by-step development of the 
TDP FDGF in Appendix [A}) These times, and the integration variables 6 and are defined by 

F{Tctoss) = Xq Tq - F{Tq) =t- X (15) 


e=(t-T + F{T) - x)/(2xo) i = {t-T-x- F{T) -f 2xo)/(2xo). 


We take /(£) = 1; the solution is then 



— (sin0 
Xo 


cos6*)-sin6l 

X 


dT 


(16) 



— (cos^ -I- sin^) 




xF{T) 


(cos^ 


sin^) 


dT. 


(17) 


The second integral extends to Tcross, the time at which the particle, moving inward, passes the potential “edge” at 
X = Xq. This integral represents the contributions to the field from the pre-crossing motion; the first integral gives 
the contribution for the motion subsequent to the edge passage. 


III. TEST TRAJECTORY 

To investigate the generation of radiation, we must choose particle trajectories for the Schwarzschild background, 
for our TDP, or for the other models that we will introduce below. There are several features we would like in a family 
of trajectories. First, we would like all trajectories in the family to conform to the analog of the physical constraints on 
trajectories in the Schwarzschild background. The most fundamental of these constraints is that dr*/dt —>■ — 1 as the 
particle approaches the horizon. For our trajectory family to be used with Eq. ([^ we therefore require dx/dt ^ — 1 
as X —>■ —oo. 

We also require that the trajectories belong to a family with at least two parameters. In this way, we can compare 
the radiation for a range of trajectories with one feature of the trajectories fixed. It is in this way that we will argue 
that the plunge radiation and QNR are not sensitive to the nature of the trajectory away from the edge. 

A third criterion for a family of trajectories is one of convenience. As a check on details, it is useful to compare 
the computation of 4/ done with the FDGF integrals, like those in Eq. 0. to the solution found by evolving initial 
data in the standard manner used in numerical relativity. For such evolution, we need initial data, and the most 
convenient initial data is static initial data, i.e., the fields that would exist at time zero if the particle had been sitting 
in a fixed position forever, and only started to move at time zero. Note that the motion before time zero is needed 
not only for evolution codes, but is also needed for the Green function solution, since the integrals, e.g., the second 
integral in Eq. 0, require a specification of the particle position back to the infinite past. For static initial data. 
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both /dt and /dt^ must vanish at time zero. This in turn means that for our trajectories both the velocity 
and acceleration must be zero at time zero. 

A desideratum, though not a requirement, is that the family of trajectories leads only to simple closed form analysis. 
The following two-parameter family for the radial motion, our “cubic trajectory” family, satisfies this as well as the 
three criteria above: 


F{T) = { ao+T-(r3 + r3)^/' r>0 
\ oo T < 0 . 


(18) 


The two parameters are oq, the initial x position of the particle, and r, the timescale on which the particle accelerates. 
For T > 0, both the velocity v and acceleration a (not to be confused with “a” used in the previous section as a proxy 
for F(t)) are 

V = -T^ (T^ -h a = -2Tt^ {T^ + t^) , (19) 

and both go to zero as T —>• 0. 

The cubic trajectory can be inverted in closed form to give T in terms of F. What is actually required for 
computations are the solution for F{T) ± T. In the F{T) — T case we must solve 


T+{t^ = p, 

where p is a known function of t,x,ao,T. The cubic equation that results has the closed form solution 

.3 \ 1/3 /,3 \ 1/3 


'' 16 64 V 4 




T = A + B 


In the F{T) F T case, the problem requires solution of 

T- (t^ F T-^ 


( 20 ) 


( 21 ) 


( 22 ) 


where cr is a known function of the trajectory parameters. By transposing one of the terms from the left to the right, 
and raising both sides to the third power, the terms cancel, leaving a quadratic equation which has the solution 


T = 


1 

2 



(23) 




FIG. 1: Comparison of computations with the Fourier domain Green function and with numerical evolution for radial infall. 
The scalar TDP test field 4' is given as a function of retarded time u. All values are given in units of the edge location xo- The 
particle generating the field moves according to the cubic trajectory of Eq. (181 with the parameters shown in each plot. The 
particle has been stationary at uq prior to T = 0. 


Figure shows the comparison between numerical evolution of the partial differential equation in Eq. ([^, and the 
evaluation of the FDGF, for radial infall, with f{t) = 1, starting from stationary initial data at t = 0 and following the 
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cubic trajectory of Eq. (181. In this figure, and in those that follow for the TDP, the scaling is set by xq, in a manner 


parallel to the scaling set in the Schwarzschild geometry, in geometrized (c = G = 1) units, by the mass parameter 
M. Thus, for example, a retarded time it = 50 is understood to mean u = SOxq and an acceleration a = 0.01 is 
understood to mean a = O.OI/xq. Note that 'h itself is also proportional to xq- 

The EDGE integrals were carried out using a straightforward Simpson’s rule scheme. The evolution scheme is the 
well-known second-order, Lax-Wendroff finite-difference scheme for second-order 1-1-1 linear hyperbolic equations; the 
delta function is represented by a narrow Gaussian, and the boundary conditions are taken to be no incoming waves 
from X = ±oo. The TDP field 4', computed, or “observed,” (i.e., numerically extracted) at a; = 3 is shown as a 
function of retarded time u = t — x. (The choice of a relatively small value of the extraction radius means that the 
results test all the terms in the TDGE, not only the terms in the radiation zone.) Although the results shown here 
are limited to the spacetime region x > Xq, and x > F(T), computations in other regions showed the same excellent 
agreement. 

Eor a scalar charged particle falling with angular motion (j) = ^(t), the source term on the right of Eq. (|^ has an 
additional factor of 5{(l) — (f>(T)). (Though angular motion in our scalar model is not representative of the change 
in the source for gravitational waves, it remains a useful check on the agreement of the EDGE and evolution.) To 
compute '^{t,x,(j)) for the i = m = 1 multipole moment we must insert /(<) = exp {±i{(j) — 4>{T)) into Eq. (11). To 
check that the EDGE integrals and numerical evolution agree for orbital motion, we take </> = 0, and we extract the 
real part of 'k, i.e., we insert cos in the integral. In order to have the angular, as well as radial motion start 
with zero acceleration, we take the angular motion to be 


Kt) = (t>o 


(r3 + r3)^/^-r 


(24) 


With this choice, orbital motion starts at (ji = 0 with zero angular velocity and angular acceleration. 
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FIG. 2: Comparison of computations with the Four ier domain Green function and with numerical evolution for orbiting infall. 
The radial motions are those described by Eq. (181, with the trajectory parameters indicated in each plot; values ar e gi ven in 
units of the edge location xo- The source has been modified with the angular factor cos (j>{T), as described in Eq. (24l. The 
examples give comparisons of three orbits that are qualitatively very different. 


Figure shows the comparison between numerical evolution and evaluation of the EDGE for orbiting infall starting 
from this stationary initial data at t = 0. The x position of the particle follows the cubic trajectory of Eq. (18) and 
the azimuthal angle is given by Eq. (24). Note that this orbital motion is not physical. The angular motion is not 
compatible with the conservation of angular momentum in the Schwarzschild background which our TDP is meant 
to model. This, however, is not relevant to the current purpose of testing the consistency of evolution and EDGE 
computations. 


IV. QNR FROM RADIAL INFALL 
A. The mathematics of QNR excitation 

To obtain insights into the plunge radiation and the excitation of QNR in the TDP, we look at a typical scenario in 
which a particle starts outside the edge of the TDP (in terms of Eq. ( [I^ , oq > a;o), crosses the edge (at a trajectory 
time Tcross)) and continues to fall in. The EDGE integrals in Eq. 0 describe this scenario. We can make this simple 




















integral even simpler by considering only the radiation zone, i.e., by ignoring the part of 'h that falls as Xjx. This 
leaves us with 




cos 6*) dT 



cos 5 


sin^ 


dT. 


(25) 


We will call the integral involving 0’s the “late” (after edge crossing) integral, and the one involving ^’s the “early” 
(prior to edge crossing) integral and we start our analysis with the late integral. The expression for 0 is illustrated 
in Fig. of Appendix in which a complete description is given of the FDGF. That figure, and the discussion 
connected to it, show that the late integral expresses the influence of the particle, consistent with causality, on the 
spacetime point at retarded time u. What is of special importance is the “damping,” i.e., the falloff of influence 
dictated by the factor e“®. 

To understand the role of this damping factor, we note that for any value of m, the damping factor is unity when 
0 = 0, that is, the contribution is maximum from the point on the trajectory labeled Tg in Appendix]^ The time Tg, 
by causality, is the latest time on the trajectory that can influence the field at retarded time u. Contributions from 
earlier points on the trajectory, at time Tg — AT, are damped by with 


A0 = (AT + [^(rg - AT) - F(Tg)]) /2xo « (1 + ug) AT/2X0 , 


(26) 


where ug = \dF/dT\, the speed of the trajectory at Tg, is a positive quantity. 

At the most basic level, this result shows that for the contribution to the late integral at any value of u, the trajectory 
is relevant only for a trajectory time less than or roughly comparable to the damping time 2xo. It also shows that 
the only property of the trajectory that matters is the velocity at Tg, unless the velocity changes substantially over a 
time comparable to one damping time, i.e., only if \d‘^F/dt‘^\ is comparable to 2a;o|fiT/(it|. 

A useful insight comes from changing the integration variable in the first integral of Eq. (25) to 0, so that the 
integral becomes 


^'late = Xq 


“““ -e (sin d — cos 0) 

6 ^ ciu , 


1 + v{0) 


(27) 


where v{9) is \dF/dT\ expressed as a function of 0 rather than T. If the velocity changes very little over the range of 
T during which most of the contribution to the integral occurs, then we can approximate v{9) = Vcmss in Eq. (27) to 
get 


^late — 


Xo 


1 T Ucross Jo 

This means that the peak of occurs at 0cross = ’’"/d, or at 


e (sin 0 —cos 0) d9 = -xq- 


‘ sin 0c 


1 + Uc 


^peak n-cross ^peak -Across T Xq 7r/4 


and that the value of d)* at the peak is 


delate,peak — 


Xo 




TT (1 


s)' 


(28) 


(29) 


(30) 


We noted that only the velocity at crossing can matter. These results show that, in fact, for the late integral the 
location of the peak of the radiation is independent of that velocity, and the strength of the radiation is very insensitive 
to the velocity. A better statement, then, is that the radiation described by the late integral is approximately 
independent of any feature of the trajectory. 

Though it is a reasonable approximation to take the velocity to be constant in the late integral, it is not useful to 
take F(T) to be a constant (in particular, to be Xg) in the second integral in Eq. (25). Graphical insights, however, 
do follow from Figure The figure shows the two sources of the radiation at the “observation” event at retarded 
time u. The 45°, “speed of light,” line from Tg to the observation event illustrates that Tg is the latest time on the 
world line that can influence that event. This, then is the graphical explanation of why the upper limit in the early 
integral in Eq. (25) is at Tg. The lower limit in this early integral is at Tcross because all of the integral in Eq. (25) 
arises from the TDGF for F(T) < xg; for T < Tcross the particle has not yet crossed that edge. 
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It is intuitively obvious that the edge at Xq must be involved in QNR. This raises a natural question about the early 
integral in Eq. (25). Since this integral is over only T < Tcross, its radiation features would seem not to be influenced 


by the existence of the edge. How can it, therefore, contribute to QNR? The question is answered graphically in 
Fig. [3 the early integral involves radiation that has interacted with the edge via the dashed path in Fig. [3 We can 
see this in the form of ^ in Eq. (@. Notice that the total distance traversed along the dashed line in Fig. |3 is the 
sum of the inward distance F(T) — Xq and the outward distance x — xq. The total time is t — T. It follows that f is 
simply the total time minus the distance traversed, so that causality requires that ^ be greater than 0. But ^ = 0 is 
equivalent to the T = Tcross upper limit of the the early integral. 


/ / 



FIG. 3: A spacetime diagram showing the influence of the particle on the “observation” point at retarded time u = t — x. 
The dashed 45° line indicates the direct, speed-of-light causality connecting the particle position at Tq with the observation 
point, and showing that Tg is the latest point on the particle world line that can influence observation at retarded time u. The 
“slower than light speed” dashed lines indicate influence from earlier times via reflection from the edge at xo- 


It is important to understand that the “slower than light” propagation involved in the early integral is not a 
manifestation of strong spacetime effects. Rather, it is associated with the fact that we are working with a multipole 
decomposition. The influence of a particle’s multipole moments on the radiation signal at an event, comes from all 
angles. The radial characteristics, then, are only an indication of the limiting speed with which multipolar information 
propagates. 


B. Models for cubic trajectories and TDP 


As the first test of the QNR excitation principle explained above, we look at a scalar charge following a radial 
trajectory in the TDP background. We use the radial trajectory function denoted by F{T) and defined in Eq. (18). 
From that equation and 


F’(ilcross) — ^0 


we have that 


Tr, 


= [(ao Tt - a;o)^ - r^] 


31 1/3 


Once Tcross is known, the velocity and acceleration at crossing can be found from 


l’(Fcross) — 


Tr, 


T=T„ 


^(Tcross) — 


dT^ 


F(T) 


T=T„ 


(r3+T„oss3)2/3’ 

2TcrossT^ 

■ (r3+r„oss3)5/3' 


(31) 

(32) 

(33) 


In all our models the velocity and acceleration are nonnegative. For simplicity, we will use the absolute values of these 
quantities in our discussions below, i.e., 


Across 

Across 


k(T'cross)|, 
I ll(-?cross) I • 


(34) 

(35) 
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FIG. 4: The radiation due to a particle in radial infall in the TDP background. The horizontal axis is the retarded time after 
the edge crossing value Mcrosa = Tcrosa — Xo- The retarded time and 4* are given in units of Xo; the acceleration is given in units 
of l/a;o. The panel on the right focuses on the region of the curves with the peak QNR. 


The radiation profiles for several radial infall trajectories, for three different values of Ucross, are shown in the two 
panels of Fig. As in Figs. and the scale of length is set by xq- All three trajectories have an acceleration 
Across = O.OOOS/xq. As a rough estimate, this acceleration multiplied by the damping time (27r divided by the 
imaginary part of the TDP QN frequency l/2xo) gives approximately Av — 0.01, on the order of a tenth of the 
crossing speed of the trajectories. It should, therefore, be a reasonable approximation to consider the infall speeds to 
be constant during the plunge and excitation of QNR. (Values of acceleration much smaller than 0.0005/xo turn out 
not to be as useful in demonstrating features of the results.) 

In the left panel the initial positive peaks are not QNR. Those large peaks, prior to QNR, are analogous to the GWs 
from the pre-QNR plunge, and we shall call these the “plunge peaks.” The QNR, a much less prominent feature of 
the profile, starts around u = Ucross, and can be seen as a small negative dip. The stark contrast of this subtle feature 
with the very visible QNR in the Schwarzschild problem is ultimately due to the fact the QNR damping frequency is 
equal to the QN oscillation frequency in the case of the TDP. In the Schwarzschild case, the QNR damping frequency 
is about a quarter of the QNR oscillation frequency. 

What is particularly important to notice is that the QNR in the right panel does not agree with the predictions in 
Eqs. (29) and (30). The QNR peak does not occur reasonably close to u — Mcross = ’’"/d, nor is the peak value of /xq 
approximately -0.19 as predicted by (30). The dependence on Ucross, in fact is markedly different from that in ( [30| ; 
while the equation predicts a peak that decreases very slightly with increasing Ucross, the results show a signihcant 
increase in the peak with increasing t’cross- This cannot be attributed to the fact that the infall speed is not exactly 
constant; in fact, the numerical evaluation of only the “late” integral for the trajectories in Fig. agrees to high 
accuracy with the predictions of Eqs. (29) and (30). Rather, we learn from Fig. that the QNR is dominated by the 
“early” integral, the second integral in Eq. (25). This is somewhat surprising, and is the hrst indication of how one 
might think about the origin of what is the dominant part of GW radiation in a binary merger. 

The question remains whether the plunge and QNR are determined only by the infall speed. In flat spacetime, after 
all, the generation of radiation is associated with acceleration. There is no hint of this in the integrals in Eq. (25), 
and it is confirmed with the numerical results in Fig. for different accelerations. With an acceleration of 0.001/^, 
for a damping time 47ra;o, there should be an increase in infall speed on the order of Av ^ 0.01. The increase in 
acceleration, with Ucross fixed implies a larger speed after crossing and a smaller speed before. In Fig. we saw that 
a decrease in both the plunge peak and the QNR excitation is associated with a smaller speed (in the case of a very 
nearly constant speed). The decrease in plunge peak and QNR excitation therefore reinforces the conclusion that it 
is the particle motion prior to edge crossing that plays the dominant role in generating this radiation. The argument 
is strengthened by noting that to very rough order of magnitude the effect of the change of acceleration is of the right 


size. 
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FIG. 5: The radiation due to a particle in radial infall in the TDP background. The horizontal axis is the retarded time 
u = t — X. (See the text for an explanation of units.) The retarded time when the first signal from the potential edge arrives is 
indicated by Tcross. Also shown in the figure is the quasinormal oscillation, the only half cycle discernible at this scale for the 
velocity and acceleration at crossing. 


V. PLUNGE AND QNR FOR TRAJECTORIES IN THE SCHWARZSCHILD GEOMETRY 

We emphasize that the simple TDP models above are exploited to give insights into the way in which the trajectories 
influence the plunge radiation and QNR. Whether they give the right insights for the black hole problem is a distinct 
question, a question we take up in this section. We use second-order Lax-Wendroff, finite-difference evolution codes [5^ 
for the Teukolsky function tpi, with our flexible cubic trajectories, to investigate whether, as in the models of the 
previous section, the plunge and QNR excitation depend only on the conditions at the “edge” of the potential. For 
this purpose, we take the edge of the potential to be at r = 3M. (The physical justification for the importance of the 
LR is discussed in Sec. VI ) 

Some comments need to be made about the way in which we judge and compare the size of QNR. Two issues arise 
here: First is the fact that unlike the ringing of normal modes, QNR has an amplitude that depends on time. We 
must therefore find a way in which the magnitude we ascribe to QNR is independent of the time at which we measure 
the amplitude. The second issue is that we are only looking at the real part of ^ 4 ; might the imaginary part contain 
different information? 

A satisfactory solution to both of these issues is to “remove” the damping from the radiation. To compare results 
for models with different trajectories, we can multiply all curves to be compared by exp (0.088965 t/M), using the 
imaginary part, 0.0889/M, of the dominant (least damped) quadrupole Schwarzschild QN mode, the mode of greatest 
astrophysical relevance. This leads to QNR having a time-independent amplitude, thereby facilitating a comparison 
of the excitation of QNR with different phases. It also shows us that the real and imaginary parts of ipi have the 
same amplitude once the damping is removed, so no information is lost by plotting only the real part of 

While this technique may be a very useful tool for working with QNR in other contexts, here we will use it explicitly 
only near the end of this section, to aid in comparing QNR with different phases. For the most part, the main value 
of this technique in the current paper is to assure us that we are not missing anything. The lessons contained in the 
curves presented below demonstrate the insights of interest, without a need for a more careful extraction of QNR 
amplitude. 


A. Models for radial infall 


We start by considering radial infall, with r* = F(t), in which F is given by Eq. (18), and in which r* is the solution 
of Eq. (§ taken to be r -|- 2MIn {r/2M — I). The models that are the basis for Fig. are chosen all to have the same 
inward speed, ulr = 0.3 at the LR. (In this figure, and all those for the Schwarzschild spacetime, numerical values 
are given in terms of the geometrized (c = G = I) Schwarzschild mass M.) Different trajectories, with this same LR 
speed, are created by varying the r parameter of Eq. (18). (The value of cq is then adjusted so that ulr = 0.3.) For a 
wide range of r values. Fig. shows that there is almost no difference in the amplitude or phase of the QNR, and that 
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there are moderate differences in the peak of the radiation at the plunge transition. Though the QNR has negligible 



FIG. 6 : The plunge radiation and QNR in the real part of the Teukolsky function 7/14 for radial cubic trajectories, all with 
vlr = 0.3 at the r = 3M LR, but with different values of the r parameter for the cubic trajectories of Eq. (18 1 . The values 
of time are shifted so that all trajectories arrive at the r = 3M LR at the same time. The values of r range from 250 to 450. 
The peak around u = 190 has a monotonically decreasing height with increasing r. Note that this peak is not pure QNR, but 
rather a transition to QNR. The curves shortly after the peak show the QNR, and demonstrate remarkable agreement in the 
amplitude and phase of the QNR for the different trajectories. 


variation, the trajectories themselves vary considerably. We illustrate this in Fig. which compares r*{T) for the 
extreme cases r = 250 and r = 450. We show in Fig. the opposite of what is shown in Fig. Here all trajectories 
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FIG. 7: Trajectories are shown for the r = 250 and r = 450 cases of Fig.[^ The to rtoi se coordinate location r* of the particle 
is shown as a funtion of particle time T, according to the cubic trajectory of Eq. (I 81 . For the r = 450 trajectory, the time 
has been shifted so that the two trajectories arrive at the r = 3M LR at the same value of T. Note that the trajectories differ 
significantly both before and after the LR crossing. 


have T = 250, while the LR speed wlr varies from 0.28 to 0.32. From the form of the trajectories, we find that 
the acceleration at LR crossing has a magnitude of approximately 2T(.rnss/'r^, which gives a value of approximately 
0.0046/M for r = 250 and 0.0026/M for r = 450. 

What insights can be gained from a comparison of the TDP results and the Schwarzschild results for radial infall? 
In the Schwarzschild case the LR plays the role of the edge. This “edge” is not sharply defined as in the TDP 
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FIG. 8 : The radiation from the plunge and the QNR in the real part of the Teukolsky function 1/14 for radial cubic trajectories, 
all with T — 250, but with different values of u = \dF/dT\ at the LR (r = 3M). The retarded time for the curves has been 
shifted, so that all trajectories pass the LR at the same time. 



FIG. 9: Trajectories for ulr = 0.28, 0.32 and r = 280, the two extreme examples in Fig. The trajectories are exactly the 
same shape since in Eq. ( |18[ ) they differ only in oq, and hence in the speed when they cross the LR. In the figure, the Ulr = 0.32 
curve is shifted so that its LR crossing time agrees with that for wlr = 0.28. 


case. Another, perhaps more important difference from the TDP case is the longer damping time. These differences 
result in an important difference in the nature of the generation of the plunge radiation and QNR. In the TDP 
case this generation was dominated by the early, pre-edge, motion. In the Schwarzschild case the generation seems 
to be dominated by the late motion. Notice in particular that in Fig. the retarded time (after crossing the LR) 
and amplitude of the QNR are almost independent of the speed at the LR. This was a specific prediction for the 
TDP provided that the QNR is dominated by the particle’s motion after crossing the potential edge. This prediction 
translates very well to the Schwarzschild case since the late (post-LR) contribution should be described by a FDGF 
like that in Eq. (13), with the TDP poles replaced by the Schwarzshild quadrupole GW poles. (To be sure, there 
will be additional functions of ui, but these will play the role of constants in the evaluation of the residues at the QN 
poles.) This implies that the late QNR excitation will be governed by an integral analogous to the late integral in 
Eq. ( [ 2 ^ . The results in Fig. [^therefore tell us that the late radiation dominates the excitation of QNR. 

Another, somewhat surprising, conclusion can be drawn from these results for Schwarzschild radial infall. Figure 
shows the plunge radiation increasing with increasing speed at the LR. This is analagous to the TDP plunge radiation 
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increasing with increasing speed at the potential edge. Insights can then be taken from Fig.[^ The plunge radiation 
is larger for larger acceleration (smaller r) at a given crossing speed. This is the reverse of what is seen in the TDP 
results, for which greater acceleration means smaller plunge radiation. Just as the TDP results indicated that the 
plunge radiation is dominated by the early (pre-crossing) motion, the Schwarzschild results indicate that the plunge 
radiation is dominated by the late, post-crossing radiation. 


B. Models for orbital trajectories 


As in the radial case, it should only be the motion near the LR that is relevant to the generation of plunge radiation 
and the excitation of QNR. This suggests that when there is angular velocity lo = d(l)/dt, it is the angular velocity 
wlr at the LR that is important, as well as the radial velocity ulr. 

If we are to investigate our hypothesis that plunge radiation and QNR excitation is determined only by the con¬ 
ditions, in this case the angular velocity, at the LR, then we need at least a two-parameter family of trajectories for 
uj(r) or a;(t), and for clarity of analysis we would like to vary the orbital motion for a fixed choice of radial motion. 
As we have explained for radial motion, the models need not be limited to families of trajectories that are physically 
correct motions. In the case of angular motion, however, we do need to consider several physical limitations on the 
nature of w. First, the orbital angular velocity d(j)/dt = w is redshifted. That is, w dies off as (1 — 2M/r) as the 
particle approaches the horizon. It turns out that we must choose orbital motions that respect this condition. If we 
do not, then the plunge radiation and QNR, for interesting values of wrr, are hidden by the much larger radiation 
coming directly from the particle motion. 

A second feature of physically correct angular motion is that duj/dt = 0 at the LR. If we were to retain this feature 
it would hamper the demonstration that it is only the value of lo at and near the LR that is important. Though this 
feature hampers our demonstration, it supports our hypothesis. Since lo near the LR is not very different from at the 
LR, our hypothesis would be wrong only if conditions on w far from the LR would be important. 

A third, and obvious, physical feature of orbits is that they must not exceed light velocity. In our notation, with 
ulr = —dr*/dt at the LR, this condition is 27Lof^^M^ < 1 — . Though this is a physical constraint on wrr it need 

not be a constraint on our exploration of the nature of the plunge radiation and QNR. 

Our choice is to take the angular motion to be 


LO = WLR 


27M2 
(1 -I- crjr^ 




(36) 


in which wlr appears explicitly as one of the parameters. For all values of the two parameters, lo in this model has 
the physically correct redshift. The parameter a governs dio/dt at the LR. For cr > 0, in our models, the angular 
frequency lo increases as the particle moves inward past the LR. 

We start the investigation of the effects of orbital motion by asking what is the interesting range of values of wrr. 
Results of evolution codes are shown in Figs. [T0| and [g For the models shown the radial speed as the LR is crossed 
is ulr = 0.3, so that the speed-of-light limit on wlr is 0.1836/M. And for all models shown in Figs. [l0| and 11 the 
a parameter in Eq. (36) is taken to be tr = 1. The results in Fig. are for very small values of wrr, values that 
do not exceed the light speed limit. For all these models, the post-plunge radiation (waveforms after t/M is 548.9) 
is QNR. For the smallest value wlr = 0.005/M the result is indistinguishable from pure radial infall. This is to be 
expected since this value is much smaller than the imaginary (damping) part of the QNR frequency, 0.088965/M. 
From the FDGF point of view this means that the motion of the source due to angular velocity is insignificant until 
the damping factor has reduced the QNR excitation integral to a negligible size. 

As wlr increases to a value comparable to the imaginary part (0.0889/M) of the QNR frequency, the orbital motion 
starts to have an effect, increasing both the plunge radiation and the excitation of QNR. For wlr = 0.01/M both can 
be seen to be distinguishable from that for wlr = 0.005/M. With a further increase of the value of cjlr, to a value 
still within the speed-of-light limitation, the effect of orbital motion becomes dramatic. 

Figure g shows the post-plunge radiation for values of wrr approaching and exceeding the light speed limit. In 
the models shown, the post-plunge radiation for both wlr = 0.05/M and 0.10/M have a period of the same order 
as the period 16.81 M of the least damped Schwarzschild QN mode. (This period is shown as a horizontal bar.) For 
wlr = 0.5/M, the very large oscillations are at a clearly different and noncostant frequency. As models progress from 
wlr = 0.1/M to higher frequency, the late time results are no longer pure QNR, but are QNR increasingly mixed 
with some higher frequency radiation. Below, in subsection |V C[ we demonstrate that for curr well above the QNR 
frequency what we are seeing is radiation coming directly from the motion of the particle. 

Orbital frequencies well above that of QNR are unphysical, so we presently limit our considerations to wrr below 0.10 
or so, and we return to the central question of this paper in the context of this subsection: Are plunge radiation and 
the excitation of QNR determined primarily by wlr, and not by the angular motion at other radii? To investigate this 
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FIG. 10: Plunge radiation and QNR excitation from models with angular frequency Mojlr =0.005, 0.01, 0.02, 0.04. All models 
have radial motion for ao = 38.6114 M and r = 600 M, and hence in all cases the particles pass the LR at time Tx = 348.90 M, 
with a speed ulr = 0.3. For all orbital motions, cr = 1. Note that each waveform shown is extracted at r* = 200 M, implying 
that the waveform at T ~ 550 originates at the LR crossing. 



FIG. 11: Plunge and post-plunge radiation for models with Mculr =0.05, 0.10, and 0.50. The radial fall is that of Fig. |10| 
The dark horizontal bar shows the period At = 16.815 M for a full wave of QNR. For Mcjlr = 0.50 the more rapid oscillations 
correspond approximately to the orbital frequency, not to the QN frequency. 


question, we consider a family of trajectories, all with the radial motion used throughout this subsection (t = 600, 
fLR = 0.3) and all with Mwlr = 0.1, a physically reasonable choice for which QNR is not dominated by direct 
radiation. We choose the models of Eq. (36) with ct = 0, 2, and -0.3. For these choices Fig. 12 shows that the time 
dependence of the orbital frequency differs significantly. The sensitivity to changes in the particle’s orbital velocity 
can be inferred from the data in Fig. mi where a change in Mwlr from 0.05 to 0.1 results in a sixfold increase in 
the amplitude of QNR. By that standard, these models should have quite different excitation of QNR if the particle 
angular velocity far from the LR were significant. 

Figure 13 shows the computed plunge radiation (real part of the Teukolsky function) for the three models. Com¬ 
parison of the amplitude of QNR cannot be made directly from plots of ip due to the phase differences in the radiation 
for the three models. We therefore apply the technique discussed at the beginning of this section, the technique of 
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“removing” the damping by multiplying by exp (0.088965 t/M). The “undamped QNR,” plotted in Fig.[^ shows the 
QNR as a series of oscillations of approximately constant amplitude. The appearance of this feature adds conhdence 
that the radiation, during this epoch, is QNR, and it greatly narrows the uncertainty about where the QNR starts. 
Most important for our current concern, it also allows a direct comparison of the amplitudes for the different models. 

From these results, with damping removed, we find that the QNR amplitudes for the cr = —0.3 model is roughly 
15% greater than that of the cr = 2 or ct = 0 models. Some variation, of course, should be expected since it is not 
precisely the ui at the LR that is crucial, but rather the uj near the LR. Here “near” means within a e-damping time 
At = M/0.0889 « IIM. The variation in the models over that time makes it surprising that there is as little variation 
in QNR excitation as is computed; the principle that the LR motion is crucial is therefore strongly supported. The 
conclusion is yet stronger for physically correct models of geodesic motion, since those models have doj/dt = 0 at the 
LR. 



FIG. 12: Particle angular velocity as a function of time for models following Eq. (361 with tJuR = 0.1, but with cr = 0, 2, —0.3. 
The a values are indicated in the plot. 



FIG. 13: The plunge radiation and QNR (real part of the Teukolsky function) for the infall models of this subsection. All curves 
here correspond to culr = 0.10, but their orbital motion, as prescribed in Eq. (361, has different values of the cr parameter, the 
values cr = —0.3, 0, and 2, as marked on the plot. Each curve has been multiplied by a factor exp (0.088965 t/M) to remove 
the QN exponential decay. This greatly clarifies when QNR starts as well as facilitates the comparison of QNR with different 
phases. 
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C. Direct radiation 


Here we demonstrate that for orbital frequencies wlr well above the frequency of Schwarzschild QNR, Muj = 
0.37367, the radiation is dominated by a field directly linked to the motion of the particle. With this we can infer 
that for frequencies near and somewhat larger than Mlo = 0.37367, the radiation will be a mixture of such direct 
radiation and QNR. We choose the specific case Mwlr = 4.0, an order of magnitude above the frequency of QNR, a 
frequency at which we should see the direct radiation completely dominate the QNR. 

To analyze the radiation at this high frequency we first need to deal with the fact that the orbital motion is so fast 
that many angular modes of radiation are generated and we lose the simplicity of a single multipole moment. For 


that reason, we present in Fig. 14 only the i,m = 2,2 part of the radiation. 

Several features of Fig. 14 are immediately interesting. The peak of the emission is around t/M « 547, the time 
at which the retarded location of the particle orbit is at the r = 3M location of the LR. The “edge” of the curvature 
potential in the Schwarzschild spacetime is, of course, not a sharp well defined edge as in the TDP, so the LR is 
only a nominal representation of an edge location. A second point to notice is that the increasing amplitude for 
t/M < 547 can be ascribed to the increasing particle acceleration, and the decreasing amplitude for t/M > 547, can 
be understood as a redshift effect. 

To delve more deeply into the nature of the “direct” radiation in Fig. we plot, in Fig. a comparison of the 
frequency exhibited by the computed wave (inferred from the spacing between the zeros of the waveform) and the 
frequency that one expects for radiation coming directly from the particle’s orbital motion. This expected frequency is 
computed by noting, for a particular “emission time” in the particle’s trajectory, the angular frequency of the particle 
as given by Eq. (36). The “observation time” is then computed for the radiation from this event to reach the wave 
extraction radius r* = 200M; it is this observation time that is the horizontal axis in Fig. The frequency is then 
multiplied by two to account for the fact that the gravitational waves are m = 2. Lastly, the Doppler shift of the 
outgoing radiation is taken into account by dividing the frequency by 1 + w, where v is the ingoing radial speed of the 
particle at the emission time. The excellent agreement between this prediction and the waveform observed leaves no 
doubt that this is “direct” radiation, and unrelated to the relatively low quasinormal frequency (indicated in the plot 
by a short horizontal line). 

The analysis of direct radiation may turn out to be more than a useful check on consistency; it may turn out 
to be an interesting diagnostic in the Kerr spacetime. In the Schwarzschild spacetime, the amplitude vs. frequency 
result has, embedded in it, information about the way in which the curvature potential acts as a high pass filter for 
radiation passing outward. It has been suggested[33] that amplitude vs. frequency curves may lead to an effective 
Kerr curvature potential. 



FIG. 14: The real part of the Teukolsky function for our infall model with iulr = 4. The l,m = 2,2 part of the solution is 
presented. 
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FIG. 15: A comparison of the frequency seen in the evolved waveform (circles) and the frequency inferred from the particle 
motion (solid curve). The waveform results above t ~ 570 are noisy due to the small amplitude of the computed waves. 
Five-point averaging was applied to smooth the waveform results below t ~ 515. The horizontal bar shows the value of a 
period corresponding to the frequency uM = 0.37367 for the least damped QN quadrupole mode. See text for details of the 
computation of the frequency inferred from the particle motion. 


VI. DISCUSSION 

Though great attention has been given to the connection of QN frequencies and the LRpfll]. much less has been 
given to the question of why the LR should play such a crucial role in the generation of plunge radiation and the 
excitation of QNR? From the point of view of multipole decomposition and the curvature potential, this is because 
the curvature potential takes a simple flat spacetime multipole form for r ^ 3M, and is negligibly small far on the 
other side of the LR (more precisely for 1 — 2M/r ^ 1). But this observation is not an explanation. An explanation 
must be rooted in the physical process being described before multipoles are projected. 

That physical process is the generation of disturbances by the motion of a particle. Disturbances propagate along 
characteristics, the same world lines as the paths of photons. Consider a particle and an observer both far outside 
(r ^ 3M) a black hole. Disturbances (scalar or whatever) of the particle will travel on the usual characteristics 
connecting the moving particle to the observer. Due to the presence of the black hole, however, there will be extra, 
non-flat-spacetime characteristics connecting the particle and observer. These will correspond to zoom-whirl orbits 
that go around the black hole, just outside the LR, and end up at the particle. Though interesting as a matter of 
principle, these extra “nonstandard” characteristics will play little role if the black hole is far away; the cross section 
for such characteristics will be too small. 

When the particle is close to the LR, however, the effect of the zoom-whirl characteristics will be dramatic. These 
“extra” characteristics will dominate over the standard characteristics as the particle passes through the LR. It is 
expected, and confirmed by studies of trajectories, that as the particle proceeds inward, its characteristics more and 
more lose the ability to reach the external observer and the particle and observer become more and more causally 
disconnected. This picture of the zoom-whirl characteristics give a fairly persuasive heuristic argument that radiation 
that will be characteristic of the black hole, plunge radiation and QNR, must come from the region of the LR. 


VII. CONCLUSION 

In this paper we have used the Fourier domain Green function to understand why the excitation of plunge radiation 
and quasinormal ringing might be determined by source conditions at the light ring of a black hole. We then verified 
that hypothesis of light ring cruciality with models both using Green function computations and numerical evolution. 
A particular insight informed by both methods was that the plunge radiation and quasinormal ringing, at least in the 
case of radial infall, can be primarily ascribed to the particle motion after crossing the light ring. 

Such insights can be viewed as a successful hrst step in understanding the plunge/ringing phase of binary black hole 
inspiral/merger, but only the first step. In particular, this work has relied on the particle perturbation linearization 
of general relativity. There are good reasons to believe that the lessons from that simplification apply to the binary 
motion of comparable mass holes, but this needs to be investigated. 

It will be necessary to extend our study to rotating holes, an extension which will be particularly interesting since 
there are two distinct light rings in the equatorial plane, one for prograde photon orbits and one for retrograde orbits, 
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and there are light rings for a continuum of radii if nonequatorial orbits are considered. Work on this analysis is 
already underway. 
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Appendix A: Fourier Domain Green Function for the TDP 


1. Method 


We give here the details of the FDGF for the simple TDP model. The method is to find a solution of Eq. (10) 
for the potential in Eq. (12). The model is simple because the homogeneous solutions of Eq. (101 are the elementary 
functions 


f e for x < Xq 

I (1 _(_ il[u}x)) , (1 — i/{ujx)) for x > xq 


(Al) 


A solution for Q is constructed by imposing outgoing boundary conditions (Q oc (1 -I- i/wx)) for x > max(a,a;o), 
and ingoing boundary conditions (G oc for x < min(a,xo). The values of x,xo and a then separate the solution 

into three forms with four coefficients. The coefficients are determined by requiring that the solution be continuous 
everywhere and have a unit discontinuity in its derivative at x = a. 

The form of the solution depends on the order of x, Xg and a, and the presentation of the solution is separated 
accordingly. We also present the results for using these FDGFs, for radial infall with f(t) = 1. 


2. Solutions for x > a >xo- 


In this case the FDGF 


^^uJ(x-a)^- 2 ^uJX 0 (^ 2 uj^xo^ + 2iujXo “ l) + e2““(l - iauj)) / i 

2aw^ (2a;2xo^ -I- 2iujxg — 1) \ xuj 


A convenient rearrangement is 


g = — 


1 


1 


1 


2ioj 2auj^ 2xuj^ 2axuj^ 


^—iLo{—x—a-\-2xQ) 


2xq{uj — uji){uj — UJ 2 ) \2iuj 2auj‘^ 2x10^ 


1 


1 1 
+ 


I 

2axuj^ 


Now we define the following two special particle times Ti and T 2 as 

Ti = t — X + F(Ti) T 2 =t — X + 2xg — F{T 2 ) . 

The meaning of these two particle times is illustrated in the spacetime cartoon in Fig. |16[ 
With these dehnitions, and with Eq. 0, we get the following: 




. ,1 1\, ^ , (t — T — X + a)^ 


dT 


(A2) 


(A3) 
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FIG. 16: A spacetime diagram showing the meaning of Ti and T 2 . The bold line represents the trajectory of the particle 
X = F{T); the dashed lines are radial characteristics. 


1 

2 




— (t — T — X — a) 

X 


{t — T — X — aY 
2ax 


dT 


rT2 




1 1 \ 2x^ 

— (cosf + sinf) + 2xo ( — I — I cosF-(cosf — sin£) 

'ax) ax 


dT. 


Here ^ is 

^ = {t — T — X — a + 2xo)f2xo . 

The first two integrals in Eq. (A41 can be combined so that the result is in manifestly finite form, 

rTi r 


(A4) 


(A5) 


1 


1 


1 rj. , ^ {t-T-x + af 

___ (,_T-x + a)- — - 


dT 



— (cos^ + sin^) + 2xq 



cos^ 


2x1 

ax 


(cos^ 


sin^) 


dT. 


(A 6 ) 


Note that this solution is good only for u = t — x less than or equal to the value of u corresponding to Merit = Across —a^Oj 
where Tcross is the time at which the particle reaches xq. For larger values of u the Green function for a > xq is no 
longer valid. 


3. Solution for inner region x < xo < a 


In this case the EDGE is 


so that 


g 


e“(“-"=)(l-iaa;) 
2a{u! — uji)(uj — UJ 2 ) 


(A7) 


4/ = — / / ———-^ dT duj 

2a[u; - u;i){uj - UJ 2 ) 

We now define T 3 by t + x = Ta + ^(Ts), as illustrated in Fig. and we define 

a={t — T + x — a)/2xo. 


With these we get 




a ^ 

• ® 

1 - -— 

sm a H-cos a 


2 xoy 

2 xo 


dT. 


a 


(A 8 ) 
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FIG. 17: A spacetime diagram showing the meaning of Ta. The bold line represents the trajectory of the particle x = F{T); 
the dashed line is an ingoing radial characteristic. 


4. Solution for intermediate region xo < x < a 


Here the FDGF is 

G = 

We define T 4 and T 5 by 




(“-“)(l-iaa;) / i \ 

2auj"^ \ (jjx) 


T 4 + F(T4) — t X 


and illustrate them in Fig. [Tsl 


4:au!'^XQ{ijj — — UJ2) V 

Ts + F{T 5 ) =t — X + 2 xo ■ 



(A9) 


(AlO) 


FIG. 18: A spacetime diagram showing the meaning of T 4 and T 5 . The bold line represents the trajectory of the particle 
X = F{T)-, the dashed lines are radial characteristics. 


For evaluation of the integrals leading to it is convenient to define 


K 4 = t — T + X — a ATs = t — T — X — a + 2xo s = Kzj2xo . 
We now break up the integral for tk into three parts and write 'F = Ii + J 2 + /a, where 




rTi r 


'■TC 


h = 


—a 


+ ^1 + —^ {—K^ + 2 xq) — — — AxqK^ + dip) 





, . , 

cose H - ^ (cose — sine) 
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n 

cose + sine — 2a;o - 
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/ — 00 

[ \a 

xy 

ax 


dT 
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(All) 


(A12) 
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(A14) 
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We can combine Ii and I 2 to give an answer that is manifestly finite 





dT 



cos £ + sin £ — 2xo 



Sig , 

— cos e H-(cos e — sin e) 

X) ax 


dT. 


(A15) 


5. Static solution for a > xq 


For comparison with evolution codes, it is important to know the initial data for evolution. This follows from 
setting w = 0 in Eq. (10) and using the potential in Eq. (12). The static solution must have the form 


= 


a/x for X > a 

PIx + ^x^ for xq < X < a 
S for X < Xq 


(A16) 


When the coefficients are matched at x = a and x = xq, so that the solutions are continuous everywhere and have a 
unit jump in their derivative at x = a, the coefficients are found to be 

a = —a^/3 — 2xg/3a ,5 = —2xg/3a 7 = —l/3a i5 = —Xg/a. (A.17) 


6. Solutions for a < xq < x 


Here the EDGE is 


2 (w - - OJ2) V 



We define Tg and 9 by 


Eg - F[Tq) = t- X 

and illustrate them in Fig. 


9 = (t — T + a — x) /2xo , 



(A18) 


(A19) 


FIG. 19: A spacetime diagram showing the meaning of Eg. The bold line represents the trajectory x = F{T) of the particle; 
the dashed line represents a radial characteristic. 


If the particle has always been moving in the a < xq region, then the Green function gives the solution 

rTe 


2 




1 2 

— (sin 9 — cos 9) -sin 9 

Xg X 


dT. 


(A20) 


In the case that mainly interests us, the particle will start at a position greater than xg and cross xg at some time 
Across- We will need to calculate 'F for values of retarded time u both less than and greater than Merit = Across ~ a;g. 
For u < Merit we find tk from Eq. (A61. For m > Merit we need Eq. (A20), but the lower limit of integration must be 
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To this integral must be 


changed from —oo, to Tcross, since the a < xq Green functio n on ly applies for T > 
added the contributions for T < Tcros s- T his comes from Eq. (A 6 | adapted to u > Ucrit- 

To understand how to adapt Eq. (A 6 |, note that when u = Ucrit both Ti and T 2 are equal to 
u > Ucrit, the case we are interested in here, both Ti and T 2 are greater than Ter, 
range of T must not exceed T = Tcross, since here the Green function is only applicable when T(T) > xq- T hus, both 
Ti and T 2 must be replaced by Tcross in Eq. (A 6 ). With this change, the first integral is missing from Eq. (A 6 I and 


Thus for 
But the integration over this 


the second integration extends only to Tcross- The complete expression then is 
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Xo 


rTe 




1 2 

— (sin 6 — cos 9) -sin 9 

Xq X 


dT 


— (cosy + siny) + 2xo ( — I — I cosy — 


2 t^ 

—(cosy — siny) 
ax 


dT. 


(A21) 


In the first integral, representing the contribution after the particle has passed inward through xq, the definition 
of Tg and 9 are, as before, those given in Eq. (A19). In the second integral, representing the contribution before the 


particle has crossed, the definitions of y is the same as that given for ^ in Eq. (A51. 


7. Solution for innner region x < a 


Here the Green function is 

■ - p—iuj{x+a — 2xo) 

g = __ g-itx{x-a) _|_ _ ^-iLd(x+a-2xo) _|__ 

2oj 2ui 2xo{uj — uji){uj — UI 2 ) 

We define T 7 and Tg as follows, and illustrate them in Fig. 

Ty + F{T^) =t + x Tg — T(Tg) = t + x — 2xo - 


(A22) 


(A23) 



FIG. 20: A spacetime diagram showing the meaning of Ty and Tg. The bold line represents the trajectory of the particle 
X = F{T); the dashed lines are radial characteristics. 


With if defined as 


ip = (t — T + x + a — 2xo)/2.xo , 


(A24) 


the solution for 41 becomes 


rTi 


rTa 


^ = -2 


'Ts 


dt — I e sin (j)dT . 


— 00 


(A25) 
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8. Solution for intermediate region xq > x > a 


Here the Green function is 


We now define ip and the times Tg, 
Tg - F(Tg) = t- X 
With these we arrive at 


^—iLo{a—x) ^—iuj{x-\-a—2xQ) 

g = -i _;; . 

2a; 4a;a;Q(a; — a;i)(a; — a; 2 ) 

Tio illustrated in Fig. 

Tio — F{Tio) = t + X — 2xo ip = (t — T + x + a 


2xo)/2xo■ 



(A26) 


(A27) 


FIG. 21: A spacetime diagram showing the meaning of Tg and Tig. The bold line represents the trajectory of the particle 
X = F{T); the dashed lines are radial characteristics. 


1 rTg 1 cTio 

41 = — - / dT — - 6““^ (cos (p + sin ip) dT. (A28) 

2 Jtio 2 


9. Static solution for a < xg 


Let a unit particle be stationary a,t x = a < xq. We decompose the static solution as follows: 

{ a/x ioi X > Xq 

P + jx ior a < X < Xq . (A29) 

S ioi X < a 

The conditions at x = Xg are continuity and smoothness; the conditions at x = a are continuity and a unit jump in 
d'^/dx. The solutions are easily found to be 


4- = 


—Xg/x for X > Xg 

— 2xg + X for a < X < Xg 
—2xg + a for X < a 


(A30) 
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